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While cosmic ray muons themselves are relatively easy to veto in underground detectors, their 
interactions with nuclei create more insidious backgrounds via: (i) the decays of long-lived isotopes 
produced by muon-induced spallation reactions inside the detector, (ii) spallation reactions initiated 
by fast muon-induced neutrons entering from outside the detector, and (iii) nuclear recoils initiated 
by fast muon-induced neutrons entering from outside the detector. These backgrounds, which 
are difficult to veto or shield against, are very important for solar, reactor, dark matter, and other 
underground experiments, especially as increased sensitivity is pursued. We used fluka to calculate 
the production rates and spectra of all prominent secondaries produced by cosmic ray muons, in 
particular focusing on secondary neutrons, due to their importance. Since the neutron spectrum 
is steeply falling, the total neutron production rate is sensitive just to the relatively soft neutrons, 
and not to the fast-neutron component. We show that the neutron spectrum in the range ~ 10- 
100 MeV can instead be probed by the (n, p)-induced isotope production rates 12 C(n,p) 12 B and 
16 0(n,p) 16 N in oil- and water-based detectors. The result for 12 B is in good agreement with the 
recent KamLAND measurement. Besides testing the calculation of muon secondaries, these results 
are also of practical importance, since 12 B (T]y 2 = 20.2 ms, Q = 13.4 MeV) and 16 N (T\/2 = 7.13 s, 
Q = 10.4 MeV) are among the dominant spallation backgrounds in these detectors. 

PACS numbers: 25.30.Mr, 25.40.Sc, 25.20.-x, 96.40.Tv 
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I. INTRODUCTION 



To reduce the cosmic-ray muon background, experi- 
ments to measure rare processes must be sited under- 
ground, where the muon flux is greatly attenuated, and 
surrounded by an active veto system, to tag the resid- 
ual muons. Even with these standard measures, muons 
are still responsible for significant backgrounds in under- 
ground experiments, via the secondary particles created 
by muon interactions with nuclei. If the muon interacts 
inside the detector, the secondary shower particles cre- 
ate unstable isotopes; some have long lifetimes, making 
it hard to associate them with particular muons. If the 
muon interacts outside the detector, it cannot be tagged, 
and "invisible" secondaries, especially neutrons, can pen- 
etrate the detector shielding. These neutrons can then 
initiate spallation reactions or nuclear recoils inside the 
detector. While these muon-induced backgrounds are al- 
ready given serious consideration at present, the next 
generation of underground neutrino, dark matter, and 
double-beta decay experiments will require both lower 
backgrounds and a better quantitative understanding of 
their characteristics. 
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Fundamental to understanding these backgrounds is 
the rate and spectrum of muon-induced neutrons [1] . The 
neutron spectrum is steeply falling over orders of magni- 
tude in neutron energy, but not uniformly so, indicating 
complexity in its formation. The total rate of neutron 
production depends primarily on the soft-neutron spec- 
trum, and has been well measured. In Table I we sum- 
marize the main characteristics of the muon flux under- 
ground at several relevant depths [2-9] . With increasing 
depth, the muon flux falls quickly and the muon average 
energy rises at first quickly and then much more slowly. 
The capture rates of neutrons produced by muons are 
also noted. The Gran Sasso rate was measured with the 
Borexino Counting Test Facility [10-12]. The rates at 
other depths were calculated using the scaling law given 
by Ref. [1] ; the result for Kamioka is fully consistent with 
the rate of 2940/ktonday measured by KamLAND [13]. 

What is needed now is a more quantitative understand- 
ing of the neutron spectrum at moderate and high ener- 
gies, as emphasized by Ref. [14]. In this paper, we con- 
sider the production of unstable isotopes as a new and 
direct probe of the moderate-energy neutron spectrum. 
We focus our attention here on the (n,p) reactions in oil- 
and water-based detectors, and show that their rates are 
a sensitive probe of the ~ 10-100 MeV neutron spectrum, 
which we calculate using FLUKA. The predicted rate of 
12 C(n,p) 12 B production is in very good agreement with 
the rate measured in the KamLAND experiment [13]. 



TABLE I: Depth, muon flux, muon average energy, and neu- 
tron capture rate at sea level, 500 m.w.e., and the Kamioka, 
Gran Sasso, and Sudbury underground laboratories. 
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p(n,7)d 




[m.w.e.] 


[u/m 2 h] 


[GeV] 


[events/kton day] 


Sea Level 





6.0xlO b 


4 


7.2x10° 


500 m.w.e. 


500 


610 


100 


8.0 xlO 4 


Kamioka 


2700 


9.6 


285 


3000 


Gran Sasso 


3800 


1.2 


320 


400 


SNOLab 


6000 


0.012 


350 


4.3 



These in-situ measurements are an important comple- 
ment to measurements made at accelerators, namely the 
experiment of Ref. [2] at CERN using muon beam ener- 
gies of 100 and 190 GeV. 

The ultimate goal of these studies is to characterize 
the muon-induced neutron spectrum precisely, at all en- 
ergies, as well as the yields of unstable isotopes pro- 
duced by muon secondaries. Here we make a step to- 
wards this goal by focusing on the reactions 12 C(n,p) 12 B 
and 16 0(n,p) 16 N in oil- and water-based detectors. We 
show below that at depths greater than a few hundred 
m.w.e., these are the only significant production chan- 
nels for 12 B in oil-based detectors and 16 N in water- 
based detectors, respectively. These are among the most 
significant spallation products in these two important 
types of detectors. By considering these single isotopes, 
with the same masses as the parents, we can isolate 
just the (n,p) production channel, and hence directly 
probe the muon-induced neutron spectrum. Even though 
the rates of these (n,p) reactions are well below the to- 
tal neutron production rates, they are still quite large: 
about 60/ktonday (calculated and measured) for 12 B in 
KamLAND [13], and 50/ktonday (calculated) for 16 N in 
Super-Kamiokande [15], both before cuts normally de- 
signed to suppress these and other spallation products. 
Both 12 B (T 1/2 = 20.2 ms, Q = 13.4 MeV) and 16 N 
(T 1/2 = 7.13 s, Q = 10.4 MeV) are unstable to (3~ decay, 
and their high production rates and endpoint energies 
make them significant backgrounds; the very long life- 
time of 16 N makes it especially pernicious. 

Section II describes our calculation of the production 
rate of the secondaries in muon showers. Section III offers 
a precise ab initio calculation of the 12 B production rate 
at different depths, and a direct comparison of our re- 
sults with the measured production rate at the Kamioka 
depth measured by KamLAND. Section IV offers a sim- 
ilar calculation for the production rate of 16 N in water. 
We draw our conclusions in section V. 



II. MUON PRODUCTION OF SECONDARIES 

Using the known muon flux underground, we used 
fluka [16] to calculate the production rates, ener- 
gies, and path lengths of all prominent secondaries, i.e., 
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FIG. 1: Cumulative path length dL(E)/dE of secondaries (in 
cm of path length per meter of fi track, and in secondary par- 
ticle energy bins of 1 MeV) generated by muons at 285 GeV, 
appropriate to the depth of Kamioka. 



7 rays, electrons (and positrons), neutrons, protons, and 
7T mesons. The fluka program is a Monte Carlo code 
able to simulate particle showers by propagating parti- 
cles according to standard interactions. The FLUKA code 
has been validated for its use in muon-induced showers 
in a number of studies. Most notably, Wang et al. [14] 
used fluka to reproduce experimental results of the pro- 
duction rate of neutrons by muons in liquid scintillator 
at several depths, and Kudryavtsev et al. [7] performed 
extensive studies on the energy spectrum and range of 
neutrons produced underground in cosmic-ray induced 
showers. 

Our FLUKA-based code was developed in the context of 
a study of the production rate of 11 C cosmogenic isotopes 
in oil-based detectors underground [12]. 

We simulated showers originating from muons at sev- 
eral relevant energies: 100 GeV (corresponding to the 
beam experiment of Ref. [2] , and also the average muon 
energy at a depth of 500 m.w.e.), 285 GeV (the aver- 
age energy at Kamioka), 320 GeV (the average energy at 
Gran Sasso), and 350 GeV (the average energy at Sud- 
bury). The use of the average muon energy should be 
adequate given that the cross sections for muon-induced 
processes scale nearly like the energy [1]. Only \x~ were 
simulated, though the results (except for muon capture) 
would be very similar for /i + . The target material in 
the simulation was the solvent of the liquid scintillator 
for Borcxino, trimethylbenzene (C9H12), with density 
0.88 g/cm 3 (incidentally, this makes up 20% of the sol- 
vent used in KamLAND [9]). The results should not 
vary greatly with other organic solvents, given that typ- 
ical values of the density and mass ratio between carbon 
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Neutron Energy [MeV] 

FIG. 2: Energy spectrum of neutrons (per meter of fj, track, 
and in neutron energy bins of 1 MeV) produced by muons at 
285 GeV, corresponding to Fig. 1. 



and hydrogen are close to those of trimethylbenzene. Ad- 
ditionally, due to the similar relevant properties, the re- 
sults for water should also be similar. We tracked muons 
for 100 meters, and for each of the prominent secondaries 
calculated the cumulative path of the particles as a func- 
tion of the particle energy, with a 10 GeV upper cutoff. 

As a representative example, in Fig. 1 we show re- 
sults for secondaries below 1 GeV produced by muons 
at 320 GeV. For each particle, this figure shows the cu- 
mulative path length dL(E)/dE traveled by all particles 
of that type at each 1 MeV bin of energy. The relative 
heights reflect both the particle multiplicities and also 
how much path length length they accumulate at each 
energy (and hence on the mechanisms of energy and par- 
ticle loss). The calculation includes all real secondary 
particles in the shower, including the abundant flux of 
bremsstrahlung photons from the muons. It is worth- 
while to note that the usually-defined "range" of the sec- 
ondary particles is not directly related to the cumulative 
path length as reported in Fig. 1: in fact, the trajectory 
of each secondary particle is broken in a large number of 
track segments, each one of them corresponding to the 
energy of the particle in that track segment; thus each 
secondary particle contributes to a large number of bins 
in the plot, from the initial energy down to lower energies 
as the particle gets slowed down along its track. 

As noted, the neutron secondaries are of special prac- 
tical importance, and here we focus just on them and the 
isotopes they produce by (n,p) reactions. Future studies 
which consider other produced isotopes will need to con- 
sider other secondaries too. In Fig. 2, we show the energy 
spectrum of the neutron secondaries. In this range, the 
neutron spectrum calculated here can be described as a 



power law ~ E~ 0,5 over ~ 10-100 MeV, and a power law 
~ E~ 2 over <~ 100-1000 MeV. Our results are consistent 
with those of Ref. [7] , which are also based on a FLUKA 
calculation. 



III. 12 B PRODUCTION IN OIL 
A. Production Reaction 

While at sea level, the capture of stopped n~ on 12 C is 
the dominant means of producing 12 B, this is no longer 
true more than a few hundred meters underground, due 
to the steeply falling fraction of stopping muons. 

The rate of stopping muons as a function of depths 
was inferred from the muon flux reported in Table I 
and from the ratio of stopping to throughgoing muons 
from Ref. [17]. At the Kamioka depth the expected 
rate of stopping muons is about 365/ktonday (this is 
consistent with the Super-Kamiokande measurement of 
220/ktonday, after taking into account the detection ef- 
ficiency of 0.65 [18]). Only negative muons can undergo 
nuclear capture, and the fraction of negative muons is 
44% [18]. The fraction of negative muons undergoing 
capture on 12 C in hydrocarbons is 7.7% [19, 20]. For 
muons undergoing capture, the branching ratio in the 
channels resulting in production of a bound 12 B state is 
18.6% [20]. Thus the expected production rate of 12 B in 
KamLAND amounts to 2.3/ktonday. 

The expected rate of 12 B production by n~ capture at 
other depths was also calculated similarly, and the results 
are summarized in Table II. It is important to bear in 
mind that beyond about 500m.w.e., the muon average 
energy and hence all secondary production rates quoted 
per meter of muon track, vary only slowly with depth. 
Accordingly, the focus of this paper is the relative rates 
of different secondary interactions. 

At shallow depths, where the 12 C(/i~, i^) channel is 
important, its rate can be very large. For example, the 
12 B rate is about 11 Hz in the inner 0.680 kton of Mini- 
BooNE (at sea level), where it is a significant background 
for the supernova detection trigger [21]. 



B. Production Reaction (n,p) 

To evaluate the contribution from (n,p) reactions, wc 
used a technique originally developed to calculate the 
production rate of the 11 C isotope in organic liquid scin- 
tillators [22], and recently exploited to calculate the 
production rate of cosmogenic isotopes in xenon detec- 
tors [23]. As shown in Fig. 1, one of the key results ob- 
tained with the fluka calculation is the cumulative path 
length traveled by secondaries of each type and energy. 
Using this, the effects of reactions of the secondaries can 
be calculated easily, without modifying their transport in 
FLUKA, provided that the reactions considered are much 
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FIG. 3: Cross sections for 12 C(n,p) 12 B and 16 0(n,p) 16 N as 
a function of the neutron energy in the lab frame. Data were 
available up to about 90 MeV for 12 C and 60 MeV for 16 O, 
beyond which we assumed the cross sections remain constant. 
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FIG. 4: Relative weighting of the 12 C(n,p) 12 B production, 
considered as an integral in dlogi? (though the data points 
are evaluated in linear steps of 1 MeV) , for neutrons generated 
in showers induced by muons at 285 GeV. 



less important than the dominant particle stopping re- 
actions. For example, for ~ 10 100 MeV neutrons, the 
(n,p) cross sections considered here are ~ 10 mb, much 
smaller than the total nuclear cross sections of ~ 1 b. 

For a secondary particle of energy E, denote the iso- 
tope production cross section by cr(E), and the appro- 
priate target density by n, so that the mean free path 
is \(E) = [na(E)]~ 1 . Thus given the cumulative path 
length dL{E)/dE calculated with fluka, the expected 
number of interactions of this type at an energy E (and 
per energy range dE) is simply \dL(E)/dE\/\(E). It is 
important to emphasize that the quantity dL{E)/dE is 
not the distance traveled by a secondary of initial en- 
ergy E; in that case, the dominant stopping reactions 
would slow the secondary and reduce its interaction rate. 
Instead, dL(E)/dE is the total amount of path length 
accumulated by all secondaries of this type, while they 
were at the energy E. 

We will indicate with Rt the total expected number 
of interactions, and hence the isotope production rate, 
given in units of per muon track length. The probability 
for each secondary to have an interaction in one of the 
channels of interest and to produce the cosmogenic iso- 
topes under study here is much smaller than unity, given 
that the cross sections for the processes of interest are 
negligible with respect to the cross sections of the domi- 
nant particle stopping reactions. Therefore we can make 
the following approximation: 



Rt 



I 



dE d jp[na(E)]. 



(1) 



is not used here directly, but only as an input to the 
second step of the fluka calculation, which handles all of 
the particle stopping reactions after having generated the 
secondaries. The relative weighting of the integral is most 
conveniently displayed with the energy on a logarithmic 
scale, i.e., in terms of dlogE ~ dE/E, the shape of this 
integrand: 



Rt ~ / dlogE 



F dL{E K(F\ 
E ^LE- a{E) 



(2) 



Note that the initial secondary particle energy spectrum 



In this paper, we consider just the (n,p) reactions for 
secondary neutrons. However, for any secondary particle, 
any reaction which can be considered as a perturbation 
to the main particle stopping reactions could be treated 
very similarly. 

Returning to the particular case of 12 C(n,p) 12 B, the 
cross section was compiled from a number of refer- 
ences [24] and is shown in Fig. 3. The same figure also 
includes the cross section for the process 16 0(n,p) 16 N, 
also compiled from a number of references [25]. 

As noted, a full simulation of muon-induced showers 
was performed with fluka [16], leading to Fig. 1. The 
product of these two figures and the energy E (i.e., con- 
sidered as an integral in dlogE) is shown in Fig. 4, show- 
ing that for this reaction, the most important neutron 
energies are ~ 10-100 MeV, probing a crucial region of 
the neutron spectrum shown in Fig. 2. For each decade 
or fraction thereof in neutron energy, the relative contri- 
bution to the integral can be immediately estimated by 
the relative height of the displayed curve. Due to the 
(n,p) cross section threshold, this reaction is insensitive 
to the very numerous soft neutrons. 
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TABLE II: Production rates for 12 B in muon induced showers 
at different depths D, given both per muon track length, and 
per volume and time. The experimental number reported by 
KamLAND [13] is also noted for comparison. 



D [m.w.e.] 


500 


2700 3800 


6000 


(E„) [GeV] 


4 100 


285 320 


350 


Process 


Rate [10" 


"%um] 




(n,p) 


1.0 10.2 


26.9 31.2 


32.2 




12.6 2.3 


0.9 0.8 


0.8 


Process 


Rate [events/ktonday] 




(n,p) 


1.4x10° 1480 


61.3 8.9 


0.1 


(m _ ,"*») 


2.1 xlO 6 390 


2.3 0.3 


0.003 


Total 


2.2x10° 1870 


63.6 9.2 


0.1 


Measured 




60 





C. Other Production Reactions 

We also examined the production channels triggered 
by 7r~ interactions: 12 B can be produced either by 
7r~ capture or by tt + photoproduction, 12 C(7,7r + ) 12 B. 
The number of tt~ produced in muon-induced showers 
at the Kamioka depth is 4.4xl0~ 3 //iin. The fraction 
of stopping pions producing 12 B isotopes in carbon is 
9.7xlCT 4 [26]. The rate of 12 B production through n~ 
capture is less than 4xl0 _6 //zm and therefore negligi- 
ble with respect to the two main channels. Concerning 
the 12 C(7,7r + ) 12 B exchange reaction, the cross section is 
<~ I fib above a threshold of 155 MeV [27], and thus the 
yield through this channel is negligible. 

Production of 12 B in organic liquid scintillators can 
also happen by interaction on the target 13 C. The low 
natural isotopic abundance of 13 C (1.1% [28]) and the 
rate of the cosmogenic reaction 12 C— > n C [22], also result- 
ing in the net loss of a nucleon from the original isotope, 
suggest that the production rates through these channels 
are negligible. 

D. Total Rates for 12 B Production 

In Table II, we summarize the 12 B production rates 
at different depths. For the four underground depths, 
the (n,p) results were obtained by direct calculations 
with fluka, as described. At sea level, the result was 
estimated by scaling the neutron production cross sec- 
tion as a ex. E Q [1]. The value chosen for the a is the 
average of the measured values on a number of unsta- 
ble isotopes produced on 12 C in the beam experiment 
at CERN: a = 0.73 [2]. The depth dependence of the 
(yu~,i^) results was obtained using the stopping muon 
fractions given by Ref. [17]. Also in Table II, the rates 
per volume Ry were obtained by 

Rv = R T <i>A M D/p)f3, (3) 

where Rt is the rate per muon track length, $ M is the 
muon flux, Mjj the detector mass, p the mass density, 



TABLE III: Production rates for 16 N in muon induced showers 
at different depths D, given both per muon track length, and 
per volume and time. 



D [m.w.e.] 


500 2700 3800 


6000 


(E») [GeV] 


4 100 285 320 


350 


Process 


Rate [lO-'V^m] 




(n,p) 


0.8 9.1 23.0 25.6 


26.3 


(m~>m) 


19.7 3.6 1.4 1.3 


1.3 


Process 


Rate [events/ktonday] 




(n,p) 


l.lxl0 D 1320 52.4 7.3 


0.07 


(m",^v) 


2.8xl0 6 530 3.2 0.4 


0.004 


Total 


2.9x10" 1850 55.5 7.7 


0.08 



and (3 a correction factor to compensate for averaging 
over the muon spectrum [2]: 

/?= 7^-4- = 0.87 ± 0.03 . (4) 

Because of how we have defined our inputs, the factor /3 
is only needed for the (n,p) calculations. 

From the values listed in Table II, one can see that the 
dominant process for the production of 12 B at depths 
greater than a few hundred meters is the (n,p) exchange 
reaction, the ((JT,^) reaction becoming much less im- 
portant. The systematic error on the production rate 
quoted in Table II due the uncertainty on the (n,p) cross 
sections is estimated to be about 5%. Other uncertainties 
and approximations probably increase this, but neverthe- 
less, the calculation is in excellent agreement with the 
rate measured at 2700 m.w.e. depth in KamLAND [13]. 
This agreement is an important confirmation of entire 
procedure for calculating both the secondaries produced 
by muons, as well as the interactions of those secondaries. 

The sea-level calculations should be taken only as 
crude estimates, since one would have to properly take 
into account the shielding of the detectors, non-vertical 
muons, unattenuated hadronic cosmic rays, etc. Addi- 
tionally, most detectors on the surface are very small, so 
that the showers induced by muons would not be fully 
contained; conversely, the small size means little shield- 
ing from interactions outside the detector. As an example 
of the importance of the 12 B production rates and mech- 
anisms, we note that the proposal of Ref. [29] to measure 
reactor i> e + e~ — > D e + e~ scattering as a test of sin 2 9w', 
the signal is a single scattered electron, and there is a 
background from 12 B beta decays [29]. 



IV. N PRODUCTION IN WATER 

Following the same procedure as above, we also calcu- 
lated the production rates of 16 N in a water-based de- 
tector. The two production channels taken into consid- 
eration are 16 0(n,p) 16 N and 16 0(^~, ^) 16 N. The cross 
section data for 16 0(n,p) 16 N are shown in Fig. 3. Note 
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FIG. 5: Yields for neutron capture, 12 C(n,p) 12 B, and 



2 C(u~ 



J B, as a function of depth, in units of per muon 



track length. In the shaded region below 500 m.w.e., the val- 
ues and their variation with depth should be taken only as 
crude estimates. 



that by comparing the cross sections on 12 C and 1 O, a 
somewhat lower range of neutron energies is relevant in 
the latter case. The fraction of stopping negative muons 
undergoing capture on 16 in water 18.4% [19], and the 
fraction of these ending in the ground state of 16 N is 
10.7% [20]. Results for the production rates of 16 N are 
given in Table III. 



V. CONCLUDING REMARKS 

In this paper we present a study of the production 
mechanism of the 12 B isotope in oil-based (organic liq- 
uid scintillator) detectors, and also of the 16 N isotope 
in water-based detectors. At depths more than a few 
hundred m.w.e. underground, their production is almost 
completely via (n,p) reactions initiated by fast muon- 
induced neutrons. We performed an ab initio calculation 
of the production rates and compared the calculated to- 
tal production rate for 12 B with data measured in Kam- 
LAND, obtaining excellent agreement. The paper offers 
a further validation of the technique exploited for the 
calculation of the rate of production of cosmogenic iso- 
topes, which was previously developed in the context of 
the study of the n C rate in oil detectors [12] and of sev- 
eral cosmogenic isotopes in xenon [23] . 

In Fig. 5 we show the variation with depth of three 
reaction rates in oil-based detectors: 

1. The 12 C(/i~, ^) 12 B rate, a measure of the muon 
flux, which is well measured and understood [17]. 



2. The p(n,"f)d rate, a measure of the muon-induced 
soft-neutron flux, which is reasonably well mea- 
sured and understood [30]. 

3. The 12 C(n,p) 12 B rate, a measure of the muon- 
induced moderate-energy neutron flux, which is un- 
certain [14]; it is quite significant that the point at 
2700 m.w.e. has been confirmed by KamLAND [13]. 

The curves for water-based detectors are similar. 

Based on these results, we note that the depth de- 
pendence of these reactions is both mild and well- 
understood. Accordingly, we place the most significance 
on the relative heights of the curves in Fig. 5. Thus 
in terms of further testing of the muon-induced back- 
grounds underground, it is difficult to make progress by 
trying to measure the mild depth dependence more pre- 
cisely. Instead, it would likely be much more fruitful to 
measure isotope production ratios at a fixed (or extrap- 
olated) depth, since these vary by orders of magnitude, 
not factors of 2. These orders of magnitude reflect both 
the strong variation of the secondary spectra with energy, 
as well as the energy dependence of the associated iso- 
tope production reactions. For example, a measurement 
of the 12 B production rate, especially relative to the to- 
tal neutron capture rate, directly probes the 10-100 MeV 
neutron flux. 

Thus, it would be very valuable if the KamLAND, 
Super-Kamiokande, and the Sudbury Neutrino Observa- 
tory experiments were to publish their detailed results 
on the relative yields of unstable isotopes produced by 
muons, as a function of distance from the muon track. 
It would be especially useful to have results on the cor- 
relations in particle yields, i.e., which isotopes (includ- 
ing neutrons) accompany each other in a given spalla- 
tion interaction, and at what distances. In the Sudbury 
Neutrino Observatory, absolute muon rates are very low, 
which restricts the possible statistics; however, since their 
intrinsic and muon background rates are so low, there is 
a unique opportunity to measure all detector activity fol- 
lowing a muon out to very large distances and times. 

The development of a well-tested physical model for 
all secondaries induced by muons would very likely allow 
more precise cuts in existing experiments, some of which 
have ~ 20% deadtime due to cuts following muons. It 
would also lead to better design considerations for future 
experiments pursuing greater sensitivity for reactor neu- 
trinos [31], low-energy solar neutrinos [32], the diffuse su- 
pernova neutrino background [34], double beta decay [33], 
and dark matter [35]. 
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